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We calculate the leading order corrections to the B — > X s ^/ decay in the standard model with 
two large flat universal extra dimensions. We find that the contributions involving the exchange of 
Kaluza-Klein modes of the physical scalar field depend logarithmically on the ultraviolet cut-off 
scale A. We emphasize that all flavor-changing neutral current transitions suffer from this problem. 
Although the ultraviolet sensitivity weakens the lower bound on the inverse compactification radius 
1/R that follows from B — > X 3 y, the constraint remains stronger than any other available direct 
measurement. After performing a careful study of the potential impact of cut-off and higher-order 
effects, we find 1/R > 650 GeV at 95% confidence level if errors are combined in quadrature. Our 
limit is at variance with the parameter region 1/72 < 600 GeV preferred by dark matter constraints. 
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I. INTRODUCTION 

The branching ratio of the inclusive radiative JJ-meson 
decay is known to provide stringent constraints on vari- 
ous non-standard physics models at the electroweak scale 
[l[, because it is accurately measured and its theoretical 
determination is rather precise. 

The present experimental world average, which in- 
cludes the latest measurements by CLEO || , Belle 3, 
and BaBar (4j, is performed by the Heavy Flavor Av- 
eraging Group [f| and reads for a photon energy cut of 
Ey > Eq with Eq = 1.6 GeV in the ii-meson rest- frame 1 

B(B -> X sl ) cxp = (3.55 ± 0.24±°;°g ± 0.03) x 10~ 4 . (1) 

Here the first error is a combined statistical and sys- 
tematic one, while the second and third are system- 
atic uncertainties due to the extrapolation from Eq = 
(1.8 — 2.0) GeV to the reference value and the subtrac- 
tion of the B — > Xdj event fraction, respectively. 

After a joint effort @, H, Q, the first theoretical es- 
timate of the total B — > X s j branching ratio at next- 
to-next-to-leading order (NNLO) in QCD has been pre- 
sented recently in Refs. [1 [nj. For E = 1.6 GeV the 
result of the improved standard model (SM) evaluation 
is given by 2 

B(B^X sl ) SM = (3.15 ±0.23) x HT 4 , (2) 



1 The very recent measurement of BaBar [6| that gives B(B — ► 
X s j) = (3.66 ± 0.85 sta t ± 0.60 syat ) x 10" 4 for E = 1.9 GeV is 
not taken into account in the average of Eq. JTJ. 

2 The small NNLO corrections related to the four-loop b — > sg mix- 
ing diagrams 3 and from quark mass effects to the electromag- 
netic dipole and current-current operator 1 1 J contributions 
are not included in Eq. p). 



where the uncertainties from hadronic power corrections 
(±5%), higher-order perturbative effects (±3%), the in- 
terpolation in the charm quark mass (±3%), and para- 
metric dependences (±3%) have been added in quadra- 
ture to obtain the total error. 

Compared with the experimental world average of 
Eq. (fTJ), the new SM prediction of Eq. |J5J| is lower by 
1.2(7. Potential beyond SM contributions should now 
be preferably constructive, while models that lead to a 
suppression of the b — > sj amplitude are more severely 
constrained than in the past, where the theoretical de- 
termination used to be above the experimental one. 

As emphasized in Refs. fl3l . [Til [Til ], among the lat- 
ter category is the model with a flat, compactified extra 
dimension where all of the SM fields are allowed to prop- 
agate in the bulk (l6l | , known as minimal universal extra 
dimensions or UED5. Since Kaluza-Klein (KK) modes 
in the UED5 model interfere destructively with the SM 
b — > S7 amplitude, the B(B — > X s j) constraint leads to 
a very powerful bound on the inverse compactification 
radius of 1/R > 600 GeV at 95% confidence level (CL) 
[l5| . This exclusion is independent from the Higgs mass 
and therefore stronger than any limit that can be derived 
from electroweak precision measurements [r7| . 

The purpose of this article is to study the phenomenol- 
ogy of B — > X~s2L m the SM with two universal extra 
dimensions [l||, El or UED6. In contrast to UED5, the 
UED6 model has additional KK particles in its spectrum. 
An interesting feature of this model is the fact that dark 
matter constraints suggest a rather small KK mass scale. 
Therefore it is very interesting to derive a bound on this 
scale from b — * 57 in UED6, taking into account the new 
KK modes. In this context, several questions will need to 
be answered: Does the leading order (LO) result depend 
on the cut-off scale, in contrast to UED5 where no cut-off 
dependence was found? If so, is this a generic feature of 
all flavor-changing neutral current (FCNC) amplitudes 
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in the UED6 model? What is the theoretical uncertainty 
stemming from the unknown ultraviolet (UV) dynamics? 

This article is organized as follows. In Sees. [Ill and 
IIIII we describe, first, the model itself and, second, the 
calculation of the one-loop matching corrections to the 
Wilson coefficients of the electro- and chromomagnetic 
dipole operators in UED6. Sec. lIVI contains a numerical 
analysis of B(B — > X s j) and the lower bound on the com- 
pactification scale 1/R in the UED6 model. Concluding 
remarks are given in Sec.[Vl In App.fA]we show how to 
compute the double sums over KK modes appearing in 
the calculation of B — > X s ^f. 

II. MODEL 

Here we briefly summarize the main features of the 
UED6 scenario. All SM fields propagate in two flat extra 
dimensions, compactified on a square with side length 
L = nR and adjacent sides being identified 20]. This 
compactification, aptly dubbed chiral square, leads to 
chiral fermion zero modes, while the higher KK modes of 
the fermions are vector-like as usual. Since the geometry 
is invariant under rotations by 180° about the center of 
the square, the model respects an additional Zi symme- 
try. It implies that the lightest KK-odd particle is stable 
and could provide a viable dark matter candidate for a 
small KK scale 1/R < 600 GeV 0. 

Solving the six-dimensional equations of motion leads 
to an orthonormal set of functions, which depend on two 
KK indices k, I corresponding to the two extra dimen- 
sions, with fc>l,Z>0orfc = Z = [18]. The model be- 
comes strongly interacting at high energy scales, so that 
it is viewed as a low-energy effective theory which is valid 
up to some cut-off scale A. From naive dimensional anal- 
ysis (NDA) (l9j |. this scale is estimated to be A ss 10/i?, 
corresponding to an upper limit Akk < k + I w 10 for 
the KK indices. 

Before electroweak symmetry breaking, all (kl) modes 
have degenerate tree-level masses mn-n = ^k 2 + P/R. 
The degeneracy is lifted by loop corrections, which lead 
to mass operators localized at the corners of the chi- 
ral square [HI, H3] ■ Additional flavor diagonal and non- 
diagonal contributions can originate from physics at the 
UV cut-off scale. Since flavor non-universal operators 
would in general lead to unacceptably large FCNC transi- 
tions, we will assume that the localized operators are fla- 



vor conserving, so that the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix remains the only source of flavor viola- 
tion. In this work, we concentrate on the leading order 
contributions from the UED6 model to B — > X s j, us- 
ing tree-level masses for those KK excitations which re- 
ceive only logarithmic corrections from loop corrections 
and boundary terms localized at the orbifold fixed points 
[19. |22|. |23|. 12-41] . This is justified since these terms are of 
one-loop order, thus leading to next-to-leading order ef- 
fects for B — ► A,,, 7. 

Upon compactification, the six-component gauge fields 
Wjfa, M = 0, . . . , 5, decompose into four-component mas- 
sive KK vector bosons W®, kl y fJ, — 0,...,3, and two 
scalar KK fields W£ 5 , kl y Here a denotes the adjoint 
group index. Following Refs. [HI, [25|, a covariant gauge 
fixing is introduced, such that W®,^ do not mix with 
^ / 45(fez)- ^ n the six-dimensional formulation, the gauge 
fixing-term reads 

£gf = - i - adiW? + d 5 W£ - g e v 6X a )] % 

(3) 

- ^7 - d'(d 4 B 4 + d 5 B 5 + g' 6 v 6X 3 )] 2 , 

where W, B are the uncompactified SU(2) and U(l) 
gauge fields with the six-dimensional gauge couplings , 
and £M are the gauge parameters. The x° are the com- 
ponents of the six-dimensional Higgs doublet 

The six-dimensional gauge couplings and vacuum expec- 
tation value are related to the four-dimensional values by 
gf = g {,) TTR and v 6 = v/R. 

The Higgs scalars mix with the fourth and fifth compo- 
nent of the gauge fields to form the would-be Goldstonc 
bosons G? kl \ of the massive vector bosons W£, kl y and 
two physical scalars a? fe « and Wjj, kl y Only the would- 
be Goldstone bosons have zero modes G^ QQ ^ , which corre- 
spond to the usual components of the SM Higgs doublet. 
For k + I > 1, the G? kl *. are dominated by the scalar ad- 
joints 5( fc ;) and B4,5(fci) while the ^//.n are composed 
mostly of the Higgs doublet elements. For the charged 
fields one finds 
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where 



X ± = 



V2 



X = W, X ,G,a,W H 



(6) 



is the tree-level squared 
The would-be Gold- 



Here M^ {kl) = m\ kl) 
mass of the W H r kl \ and 
stone bosons recerve the unphysical squared mass 

^M^y^n from gauge fixing. Similar expressions hold for 
the neutral fields, taking into account a small mixing be- 



tween Wf kl) 



and Bn.n. However, since they do not con- 
tribute to the process B — > X s j at leading order in the 
electroweak interactions we do not give them here. 

As mentioned above, the masses of the KK modes re- 
ceive corrections from loop and UV effects, which are de- 
pendent on the cut-off scale A. Since G^j,^ and W^ H ^ 
are protected by gauge invariance, the dependence on 
A is only logarithmic [22|, so that the mass corrections 
are small compared to 1/R and can be neglected in a 



LO calculation. The 



\ki) 



scalars, however, can receive 



contributions proportional to A 2 to both their bulk and 
boundary mass terms [24]]. 

In order to obtain a small mass term for the zero mode 
Higgs doublet, the bulk and boundary mass terms need 
to be tuned to cancel to a large extent. However, inde- 
pendent of this tuning, the higher KK modes can receive 
sizeable contributions from these terms. As a result, the 



scalars can be heavier or lighter than the other parti- 
cles of the same KK level. 3 We include the A 2 corrections 
to the masses based on the following parametriza- 
tion of the UV-induced mass terms: 



C D 



— (d(x 4 )S(x 5 ) + S(L - x 4 )6(L - x 5 )) m\ x 



(7) 



Although the UV physics is not specified, these mass pa- 
rameters are expected to stem from loop contributions of 
the UV dynamics, so that 



16tt 2 



rA 2 = 



N 2 
JV KK 



16tt 2 R 2 



i = l,2, 



(8) 



with hi 2 = O(l). Using the explicit form of the KK 
wave functions from Refs. [H, [2(| and tuning the bulk 
mass m 2 H bulk to exactly cancel the A 2 correction to the 
zero mode of the Higgs doublet, the masses of the 
scalars are found to be 
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FIG. 1: One-loop corrections to the b —* s'y amplitude in the 
UED6 model involving the KK modes of the would-be Gold- 
stone, G,, n , the W-boson, W and the scalar fields a,,_, 



p{kl) ' 



(fcl) 



and W^ kl y Diagrams where the SU(2) quark doublets Q( fei ) 
are replaced by the SU(2) quark singlets ti^n are not shown. 
Here i = u, c, t. See text for details. 



We will estimate the theoretical uncertainty from the un- 
specified UV physics by varying the coupling constants 
h\^2 of the boundary mass terms independently in the 
range [0, 1] which corresponds to either decoupling or 
strong coupling. 

The boundary mass terms could cause mixing among 
KK modes and one would need to re-diagonalize the mass 
matrix to find the eigenstates if they are large. To have 
a light Higgs boson, we assume that these mixing mass 
terms are tuned to be much smaller than 1/R, so that 
we can treat them as small perturbations and ignore the 
higher-order mixing effects. 

The small KK scale suggested by dark matter con- 
straints would lead to interesting signals at the Fer- 
milab Tevatron and the CERN Large Hadron Collider 
19 , [2(| [2?| as well as the International Linear Collider 
28]. However, strong bounds on the compactification ra- 
dius can arise from heavy flavor physics. In particular, 
the FCNC decay B — > X s j, which shall be studied in the 
following, is known to put stringent constraints on vari- 
ous beyond the SM physics scenarios at the electroweak 
scale. 



III. CALCULATION 



3 This problem already arises in UED5, but was not discussed in 
previous analyses of B — > X s j for this model [13|. Il4t [l5| . 



We work in an effective theory with five active quarks, 
photons and gluons obtained by integrating out the elec- 
troweak bosons, the top quark, and all the heavy KK 
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modes. Adopting the operator basis of Ref. [29(, the ef- 
fective Lagrangian relevant for the b — ► sj(g) transitions 
at a scale \i reads 



4G F ^ ° 
£ e ff =£qedxqcd H 7=rV^*Vf& 2 Cj[p)Qj 

v 2 1=1 



(10) 



where the first term is the conventional QED and QCD 
Lagrangian for the light SM particles. In the second term 
Gf and denotes the Fermi coupling constant and the 
elements of the CKM matrix, respectively, while Ci(//) 
are the Wilson coefficients of the corresponding operators 
Qi build out of the light fields. Terms proportional to the 
small V u b mixing, which will be included in our numerical 
results, have been neglected above for simplicity. The 
same refers to higher-order electroweak corrections (30j |. 

The operators Q\ & are the usual four-quark opera- 
tors whose explicit form can be found in Ref. [29|. The 
remaining two operators, characteristic for the b — > sj{g) 
transitions, are the dipole operators 



Qs = ^(sL^T a b R )G«„ . 



(11) 



Here e (g) is the electromagnetic (strong) coupling con- 
stant, qL,R are chiral quark fields, F^ v ((?„„) is the elec- 
tromagnetic (gluonic) field strength tensor, and T a are 
the color generators normalized such that Ti(T a T b ) = 
<5 a h/2. The factor m& in the definition of denotes 
the bottom quark MS mass renormalized at /i. 

The relevant quantity entering the calculation of 
B(B — ► X s ~f) is not CV(/i) but a linear combination 
Cj (/i) of this Wilson coefficient and of the coefficients of 
the four-quark operators. The so-called effective Wilson 
coefficients relevant for b — > sj(g) are [3l[ 



for i = 1, ... , 6, 



^{p) + E Vj C j(tJ') fori = 7, 



3=1 
6 



(12) 



C 8 (m)+ HzjCjin) for l = 8, 

3=1 



where yj and Zj are chosen so that the LO b — > s"f(g) ma- 
trix elements of the effective Lagrangian are proportional 
to the LO terms C^°\fi). In the MS scheme with fully 
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anticommuting 75, one has y = (0,0,—^,— 
&ndz= (0,0,1, -|, 20, -f) 

We further decompose the effective coefficients into a 
SM and a new physics part 



Cf( M ) = C^ M (M)+ACf(M), * = 



(13) 



and expand the latter contribution in powers of a s as 
follows 



OO 
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ACf ( "V). 



(14) 



In the case of UED6, new physics affects the initial 
conditions of the Wilson coefficients of the operators 
in the low-energy effective theory while it does not in- 
duce new operators besides those already present in the 
SM. To find the LO corrections from the UED6 model 
to B{B — ► Xg'y) one has to consider all the one-loop 
one-particle-irreducible diagrams contributing to the pro- 
cesses b — > S7(#). The one-loop b — > s-f diagrams are 
shown in Fig.[TJ Before performing the loop integration, 
the Feynman integrands are Taylor-expanded up to sec- 
ond order in the off-shell external momenta and to first 
order in the bottom quark mass. Thereby only terms 
which project on Q7 after the use of the equations of 
motion arc retained. The calculation for the b —t sg am- 
plitude proceeds in the same way. The relevant Feynman 
rules have been derived from Ref. 1 1811 and implemented 
into a model file for FeynArts 3 [32J, which has been 
used to generate the necessary amplitudes. At tree-level, 
the interactions between SM and KK fields preserve both 
KK numbers. Consequently, only diagrams where all par- 
ticles in the loop have the same KK index (kl) have to 
be taken into account. 

At the matching scale /io = 0(m t ) the LO results for 
the UED6 initial conditions read 



AC 



cff(O) 



(Mo) 



for i = 1, 

-lJ2'A ( °Hx kl ) forz = 7, 



,(». 



k.l 



(15) 



iE'F {0) M for* = 



k.l 



where the / superscript in the summation indicates that 
the KK sums run only over the restricted range k > 1 
and I > 0, i.e. J2' k ,l = T, k >i Ez>cr 

We decompose the Inami-Lim functions as 



E *J°W 



X = A,F, 



I=W,a,H 



(16) 



where the function xff a H (xki) describes the contribu- 
tion due to the exchange of KK modes of the would-be 
Goldstone, GvjU, and the H^-bosons, W^ k[ y the scalar 

fields af M) and W± {kl) . Here x kl = (fc 2 + 1 2 )/{R 2 M^). 

Our results for the LO Inami-Lim functions entering 
Eq. (fl6|) are given by 
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x t (-6x1, + (6x| - 9x 4 - 9) x fci + (7 - 2x t )x t - ll) 



+ -^x k i(xki + 1) In 



12(x t - l) 3 



{x k i + l)(x t + x k i((x k i - x t + 4)xi - 1)) / Xjy + x t 
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2(x t - iy 



Xkl + 1 



x t (6 (xf - 3x t + 3) x^ - 3 (3xf - 9x 4 + 2) x k i - 7x 2 + 29x t - 16) 



36(x t - l) 3 



+ a Xkl ( x kl ~ Xkl ~ 2 ) ln ( ^ 1 
6 \Xki + 1 



(x fci + 1) (xl, + (4x t - 2)x fei + xt(3x t -2)) / x k i + x t 

In ' 



6(x t - l) 4 

x t (6 (xf - 3xt + 3) x 2 ., + 3 (3x| - 9x* + 10) x k i + 2x 2 - 7x t + ll) 



xu + 1 
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2(xt - l) 4 

r 



In 



V x k i + 1 



(17) 



(18) 



(19) 



(20) 



(21) 



(22) 



Here Xf 



TO 2 (/io)/M^. Our results for the 



Xyy ■ (xki) + Xa \xm), X — A,F, agree with the expres- 
sions for the one- loop dipole functions given in Ref. [3]. 
We note that there is a misprint in the last line of 
Eq. (3.33) of the latter paper. Obviously, the term 
ln((x„ +x t )/(l + xt)) should read ln((x n +x t )/(l + x n )) 
with x n — n 2 /(R 2 M'^ v ) and n the single KK index ap- 
pearing in UED5. 

For the numerical analysis, the results in Eqs. (jTTJ) to 
(|22|) need to be summed over the KK indices k,l. This 
summation can be performed analytically employing an 
expansion for large l/R, as explained in App.[A] For 
zero boundary mass contributions, hi t 2 = 0, we obtain 
the following approximate formulas 



k.l 



0.686134+ 0.162912 Ax t 



0.316677+ 0.075190 Ax t 



(23) 
(24) 



k.l 



£V(°>(x w ) 



k.l 



E'^SW 



k.i 



k.l 



— ^\n(A 2 R 2 ) 
288 x y ' 

0.86695 + 0.205844 Ax t 



0.791563 + 0.187945 Ax t 



0.211118 + 0.050127 Ax t 



0.158339 + 0.037595 Ax t 



(25) 



(26) 
(27) 

(28) 



k.l 



where x = \ j{R 2 M 2 / ) and Ax t = x t - (165/80.4) 2 . Note 
that in the above formulas we have only kept the leading 
terms in the 1/x expansion for simplicity. The coeffi- 
cients of the logarithms in Eqs. (|25p and (|26p are exact 
in the limit of an infinite number of KK modes. We em- 
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phasize that the given approximations are for illustrative 
purpose only. In our numerical analysis we will through- 
out employ the exact double series ^2 k l Xj°\xki), X = 
A,F,I — W, a, H, summed over the restricted range 
k > 1, I > 0, and / + k < N KK . 

We see from the latter equations that while the one- 
loop Gf k[) and W± H 

(ki) corrections to ACy^ (0) (fi ) are 
insensitive to the UV cut-off scale A or, equivalently, 
-/Vkk: the contributions due to a, fe « exchange depend log- 
arithmically on A. The different convergence behavior 
is closely connected to the unitarity of the CKM matrix 
which results in a Glashow-Iliopoulos-Maiani (GIM) sup- 
pression [33| of the higher KK mode contributions to the 
double sums in Eqs. to (f2"g)) . In the case at hand, 
the GIM mechanism leads to a hierarchy of the various 
contributions to ACy^ ' (/in), with X^) H (xki) propor- 
tional to l/(fc 2 + I 2 ) 2 and X ( a\x k i),X = A,F, scaling 
like l/(fc 2 + I 2 ) for large values of I, k. The extra power 
of k 2 + I 2 in the contribution from diagrams with o.^ k ^ 
exchange, that leads to the logarithmic divergent results, 
stems from the left- (right-handed) top quark Yukawa 
coupling enhanced part of the a^ kl ^U^ kl ^b (a^^sU*^) 
tree-level vertex. No such terms are present in the flavor- 
changing vertices involving and W^ H , k n . 

The logarithmic divergences appearing in Eqs. (|25p and 
(l2"fj)) would be cancelled by counterterms at the scale 
A at which perturbativity is lost in the higher dimen- 
sional theory. Our calculation only determines the lead- 
ing logarithmic corrections associated with the renormal- 
ization group (RG) running between A and 1/R. The 
corresponding initial conditions contain incalculable fi- 
nite matching corrections from the unknown UV physics. 
Assuming that the RG effects dominate over the finite 
matching corrections and that the UV completion of 
the UED6 model has a CKM-type flavor structure, the 
UV sensitivity can be absorbed into a logarithmic de- 
pendence on AR or, equivalently, -ZVkk- To gauge the 
theoretical uncertainty associated with the unknown UV 
completion we will vary Nkk in the range [5, 15] around 
A K k = Ai? w 10 as estimated by NDA. The choice of the 
lower value of -/Vkk is motivated by the observation that 
for Nkk < 5 the non-logarithmic terms in ; X^ (xki ) , 
X = A, F, become numerically of the same size as the 
logarithmic ones. Since the choice of the upper value of 
A^kk has no impact on our conclusions we choose it sym- 
metrically. We mention that the requirement of unitarity 
of gauge boson scattering at high energies [34J generically 
leads to values of Ai? notably below the NDA estimate 
N KK a 10. 

The individual contributions A/Cy^' ' (p-o), I — 
W, a, H, to the UED6 initial conditions of the dipole op- 
erators as a function of 1/R are shown in Fig. [5] The 
contribution due to the exchange of and W^ kl ^ 

and W^( H ) (green/medium gray) KK modes are de- 
picted as yellow (light gray) and green (medium gray) 
curves, while the red (dark gray) bands and the black 
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FIG. 2: A/Cy g (/^o) as a function of 1/R. The different curves 
correspond to the individual contributions due to the exchange of 
KK modes of the would-be Goldstone, G^r ,, , and the IF-bosons, 

^M(fcz) (y e ll° w /light g ra y). the scalar fields a^ ; j (black) and 
VFff(fci) (green/medium gray), respectively. The lower (upper) 
borders of the red (dark gray) bands correspond to TVkk = 5 (15) 
while the black lines represent the results for -/Vkk = 10. See 
text for details. 



lines illustrate the a, fc « corrections. The lower (up- 
per) borders of the red (dark gray) bands correspond to 
Akk = 5 (15) while the black lines represent the results 
for A^kk = 10. We see that in both cases the contribu- 
tion involving o/n,n exchange is by far dominant and its 
variation with Nkk is non-negligible. Nevertheless, the 
large positive corrections to AC^ \p,o) already start 

to exceed the SM values Cy^gj^^o) ~ —0.19,-0.10 in 
magnitude for 1/R a 240, 335 GeV in the most conser- 
vative case Akk = 5. The observed strong enhancement 
of the initial conditions Cy^°^ (/in) will play the key role 
in our phenomenological applications. 4 



4 For compactification scales 1/R 100 GeV it would even be 
possible to reverse the sign of Cy ff (^tf,) with respect to its SM 
value C| ff (/ii,) ~ —0.37. This possibility is disfavored on general 
grounds by the experimental information onS-t X s Z+i — [35t . 
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FIG. 3: B(B -> X s >y) for E = 1.6 GeV as a function of 1/R. 
The red (dark gray) band corresponds to the UED6 result. The 
68% CL range and central value of the experimental/SM result 
is indicated by the yellow/green (light/medium gray) band un- 
derlying the straight solid line. See text for details. 



Another main observation of our work is, that in the 
UED6 model the Z- (AC), photon (AD), gluon pen- 
guin (AE), and the |AF| = 2 boxes (AS") all behave 
as l/(k 2 + I 2 ) for large values of k,l. In contrast, 
|AF| = 1 boxes (AB^^u) show an asymptotic l/(k 2 +l 2 ) 2 
behavior after GIM. The corresponding UED6 Inami- 
Lim functions therefore exhibit the following behavior: 
AC, AS" cx x 2 /x\n(A 2 R 2 ), AD, AE cx x t /x \n(A 2 R 2 ), 
and AB vl/ ^i cx Xt/x 2 . This implies that the logarith- 
mic cut-off sensitivity first seen in Eqs. (f25|) and ([26]) is 
a generic feature of all FCNC transitions in the UED6 
model. A dedicated study of neutral meson mixing, rare 
K- and B-decays in UED6 is left for further work. 



IV. NUMERICS 

The UED6 prediction of B(B -> A s7 ) for E Q = 
1.6 GeV as a function of 1/R is displayed by the red (dark 
gray) band in Fig. [3] The yellow (light gray) and green 
(medium gray) band in the same figure shows the experi- 
mental and SM result as given in Eqs. |T]) and respec- 
tively. In all three cases, the middle line is the central 
value, while the widths of the bands indicate the uncer- 
tainties that one obtains by adding errors in quadrature. 
The central value of the UED6 prediction corresponds 
to Nkk = 10 and hi 2 = 0. The strong suppression of 
B(B — > X s -f) in the UED6 model with respect to the 
SM expectation and the slow decoupling of KK modes is 
clearly seen in Fig. [31 

In our numerical analysis, matching of the UED6 Wil- 
son coefficients at the electroweak scale is complete up to 
leading logarithmic order, while terms beyond that order 
include SM contributions only. For the reference values of 
the renormalization scales [1q, /Zb, fj, c = 160, 2.5, 1.25 GeV, 



we utilize the formula 

X sl ) = [3.15 ±0.23 

- 8.03 ACf (0 Vo) - 1-92 ACf {0) (^) 

+ 4.96 (AC c 7 m M) 2 + 0.36 (ACf ( °V )) : 

+ 2.33ACf' (0) (^o)AC 8 cff(0) ( A/0 )l x HT 4 , 



(29) 



which has been derived based on the NNLO SM results 
of Refs. [1, Qil HI] . For the remaining input parameters 
we adopt the central values and error ranges that can be 
found in Ref. 0. 

The theoretical uncertainty in the UED6 model is esti- 
mated by scanning Vkk, the couplings h±^ of the bound- 
ary mass terms, and the matching scale /in in the range 
[5,15], [0,1], and [80, 320] GeV for the largest possible 
variations. The combined theory error does not exceed 
±g 7 % for 1/R in the range [0.4, 2.0] TeV. Larger rela- 
tive errors of above appear for 1/R = 300 GeV. 
Whether the quoted numbers provide a reliable estimate 
of the cut-off and higher-order corrections to B(B — > 
X s ~f) in the UED6 model can only be seen by perform- 
ing a next-to-leading order (NLO) matching calculation. 
Such a calculation seems worthwhile but is beyond the 
scope of this work. The parametric uncertainty due to 
the error on the top quark mass is below i 3 % for 1/R 
in the range [0.3, 2.0] TeV and thus notably smaller than 
the combined theory uncertainty. 

Since the experimental result is at present above the 
SM one and KK modes in the UED6 model necessarily 
interfere destructively with the SM b — ► 57 amplitude, 
the lower bound on 1/R following from B(B — > X s j) 
turns out to be much stronger than what one can derive 
from any other currently available direct measurement 
[26j . If experimental, parametric, and theory uncertain- 
ties are treated as Gaussian and combined in quadrature, 
the 95% CL bound amounts to 650 GeV. In contrast to 
the upper limit coming from the dark matter abundance 
the latter exclusion is almost independent of the Higgs 
mass because genuine electroweak effects related to Higgs 
boson exchange enter B(B — > X s j) first at the two-loop 
level. In the SM these corrections have been calculated 
[3fl ] and amount to around —1.5% in the branching ra- 
tio. They are included in Eq. (f2"9")l . Neglecting the cor- 
responding two-loop Higgs effects in the UED6 model 
calculation should therefore have practically no influence 
on the derived limits. 

The upper (lower) contour plot in Fig.Q]shows the 95% 
CL bound of 1/R as a function of the experimental (SM) 
central value and error. The current experimental world 
average and SM prediction of Eqs. (fT]) and @ are indi- 
cated by the black squares. These plots allow to monitor 
the effect of future improvements in both the measure- 
ments and the SM prediction. Of course, one should 
keep in mind that the derived bounds depend in a non- 
negligible way on the treatment of theoretical uncertain- 
ties. Furthermore, the found limits could be weakened by 
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FIG. 4: The upper/lower panel displays the 95% CL limits on 
1/7? as a function of the experimental/SM central value (hori- 
zontal axis) and total error (vertical axis). The experimental/SM 
result from Eq. (JTJ/Eq. (01 is indicated by the black square. The 
contour lines represent values that lead to the same bound in 
TeV. See text for details. 



the NLO matching corrections in the UED6 model which 
remain unknown. 



V. CONCLUSIONS 

We have calculate the leading order corrections to the 
inclusive radiative B — > X s j decay in the standard model 
with two universal extra dimensions. While the one- 
loop matching corrections associated to the exchange of 
Kaluza-Klein modes of the would-be Goldstone, Gfkiy 

the W-boson, W^, k n , and the physical scalar W^, kl s are 
insensitive to the ultraviolet physics, we find that contri- 
butions involving o^ kl ^ scalars depend logarithmically on 
the cut-off scale A. We have emphasized that in the con- 
sidered model all flavor-changing neutral current transi- 
tions suffer from this problem already at leading order. 
Moreover, we have included formally next-to-leading, but 
sizeable mass corrections to the Kaluza-Klein scalars that 
depend quadratically on the scale A. Although the ultra- 
violet sensitivity weakens the lower bound on the inverse 
compactification radius I / R that can be derived from the 
measurements of the B — ► X s j branching ratio, a strong 
constraint of 1/R > 650 GeV at 95% confidence level is 
found if errors are added in quadrature. Our bound ex- 
ceeds by far the limits that can be derived from any other 
direct measurement, and is at variance with the parame- 
ter region preferred by the dark matter abundance. This 
once again underscores the outstanding role of the inclu- 
sive radiative £?-meson decay in searches for new physics 
close to the electroweak scale. 
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APPENDIX A: EVALUATION OF KK SUMS 



Here we show how to approximate the double sum over KK levels (kl) appearing in Eq. (fTS"]) . Following Ref. [If 
we first introduce the integrals 



I n (a) = (-l) n a 



n n+1 



dy — : 
o ay 



Xkl 



(Al) 
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where n = 0, 1, . . . , and Xki = (k 2 + l 2 )x with x — 1 / (R 2 M^) ■ Obviously, I n (0) = 0. These integrals allow use to 
express the logarithms appearing in Eqs. (jTTJ) to as 



In 



xh In 



xi. In 



4z In 



xkl 


4- n 


Xkl 


+ 1 


Xkl 


+ a 


Xkl 


+ 1 


Xkl 


+ a 


Xkl 


+ 1 


Xkl 


+ a 


Xkl 


+ 1 



Jo(a)-/ (l), 

h(a) - h{l) - 1 + a , 



h{a) - J 2 (l) + ^ - Xfe; + x H a - ^a 2 . 



(A2) 



/.-)(") - /.-id • - ^ + ^arw - x 2 kl +x 2 kl a- ^x u a 2 + ^a 3 . 



with a = or x t . We note that Eq. (D.3) of Ref. [25| is missing an overall minus sign on its right-hand side. 

Since the individual building blocks I n (a) behave as l/(k 2 + I 2 ) for large fe, I, the corresponding double series over 
the KK levels diverge logarithmically. We regulate the appearing divergence analytically 



with S > 0. Then one has 



J» = (-l)"a 



oo oo „i 

J2'lfa) = (-l) n a n+1 j:j: / dy 

k,l k=l 1=0 



n n+1 



o dy (aj/ + x fci ) 1+5 ' 



(A3) 



(ay + 



1+* 



(— l)"a™ +1 °° 00 r 1 



r(i- 




(-i)" 


a™ +1 


4T(1 


+ ^) 


(-1 


P 


4T(1 - 


ho) J 



fe=l Z=0 ' 



(A4) 



1 e 



-ayt 



dy y n jf £ ftt«(i?3(0,c- !rt y 

dtt -i-n+s r §3 ^ e -xt^ _ ^ ^ + ^ _ r ^ + ^ at ^ ^ 



where in the first step we have used the Mellin-Barnes representation 

1 1 



T(1 + S)J 



dtt 5 e~ st 



(A5) 



Here <& 3 (u,q) = 1 + 2 £^ =1 g" 1 ' cos(2tou), r(z) = J^dtt^e^ and T(u,z) = f^dtt^e'*, denotes the elliptic 
theta, the Euler gamma, and the plica function, respectively. 

The integration over t in the last line of Eq. (|A4|) cannot be performed analytically. Yet using 



03 (0,e" z ) pb < 



for z < t/tt, 
l + 2^e-™ 2z , forz> % /¥, 



n+1 



(A6) 



m— 1 



and expanding the integrand in powers of \jt in the latter case, we can perform the integration piecewise and 
approximate the double series as 



£j£(o)»Z£(a) + Mo). 
k,l 



(A7) 
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The integration over t £ [0, t/tt/x] leads to the relatively compact formulas 



(-l)™7ra n+1 1 
4(n+l)x 5 



1 

8^ 



2^xE 2 ( ^— ) - x ( 2T (o, ^ ) + In 



+ < 



+ 2 (7ra(l - In a) - (VtF + 7e) x) 
(-1)" 



for n = 0, 



An{n + l) 2 x 



(n+ l)^a n - a" x(n+l) 



(A8) 



+ 7ran(n + 1) [T 0, 



+ ln(a) I — nan 



+ (n + 1) (n (1 - + 1) n-' n/2 x H+1 



x r(n + 1) 



for n = 1,2, 



where we have expanded the result around S — and dropped all terms that vanish in the limit S — > 0. Furthermore, 
E m {z) — dtt~ m e~ zt and 7e ~ 0.577216 is the exponential integral function and the Euler constant. 

The integration over t S (^/n/x, oo) is finite in the limit <5 — > 0. For all double sums y\ ; I n (a) appearing in Eq. (|A2[) 
we were able to find analytic expressions. Since the results turn out to be rather lengthy and not very informative 
we refrain from giving them here. Short numerical expressions for the h n (a) can be obtained in the large x limit. 
Keeping terms up to third order in 1/x, we find 



0. 



h n (a) 



Combining Eqs. 



k,l 



(-l)™7ra" +1 
4(n + l)x 



0.184616 a 


0.25221a 2 

2 + 


0.259202 a 3 


for n 


X 




X z 


a; 3 




0.0923082 a 2 


0.16814 a 3 


0.194402 a 4 


for n 


X 




X 2 


X 3 




0.0615388 a 3 


0.126105 a 


4 0.155522 a 5 


, for n 


X 




X 2 


X 3 




0.0461541a 4 


0.100884 a 

+ 2 


5 0.129601a 6 


, for n 


X 




X z 


X 3 




finally arrive 


at the following large x approximations 






0.644381 a 


0.751902 a 2 0.387481a 3 






X 


X z 


X 3 






0.322191a 2 


0.501268 a 3 


0.290611a 4 






H 

X 


X 2 


X 3 


- In x^j + < 














0.214794 a 3 


0.375951a 4 


0.232489 a 5 






X 


X 2 


X 3 






0.161095 a 4 


0.300761 a 5 


0.193741a 6 






4 

X 


X 2 


X 3 



(A9) 



for n = 0, 



for n = 1, 



(A10) 



for n 



for = 3. 



The term 1/5 — lux in Eq. (TATOl) implies that one should include counterterm contributions from physics at the UV 
cut-off scale A that cancel the divergences. Our calculation only determines the RG running contribution between A 
and 1/R, given initial conditions at A. Assuming that the unknown finite matching corrections are small and have a 
CKM-type flavor structure, the divergences can be absorbed into a cut-off dependence by switching from analytic to 
cut-off regularization employing the approximation 



In a; « ln(A 2 i? 2 ) , 



(All) 
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with A not much larger than 1/R. We remark that the latter assumptions are self-consistent because the finite 
matching corrections are formally of next-to-leading logarithmic order. 
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